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ABSTRACT 

The  results  of  seven  different  tests  of  a  single  stage 
axial  turbine  with  converging-diverging  stator  nozzle  are 
reported.   From  measurements  in  a  special  test  rig,  losses 
occurring  in  the  stator  and  rotor  blade  rows  were  separately 
calculated  and  the  performance  of  the  stage  was  also  deter- 
mined.  The  rotor  speeds  varied  from  9,500  to  18,600  r.p.m. 
and  the  pressure  ratios  varied  from  1.75  to  3.25.   The  Mach 
numbers  at  the  stator  exit  varied  from  0.79  to  1.38.   The 
results  for  this  turbine  are  appraised  and  a  procedure  is 
demonstrated  for  smoothing  loss  coefficient  data  from  the 
turbine  rig.   Test  rig  improvements  reported  include  the 
design  and  construction  of  a  new  flow  nozzle  and  the  revision 
of  the  data  reduction  programs  to  access  a  Hewlett-Packard 
Model  9867B  Mass  Memory  unit. 
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I.   INTRODUCTION 

The  transonic  turbine  test  rig  in  the  Turbopropulsion 
Laboratory  at  the  Naval  Postgraduate  School  was  designed 
to  determine  the  effect  of  blading  design  on  turbine  effi- 
ciency and  to  allow  the  separate  determination  of  rotor  and 
stator  losses. 

Preceding  the  work  of  Solms  [Ref .  5]  in  which  several 
different  turbine  geometries  were  tested,  the  separation  of 
the  losses  in  the  blade  rows  was  not  achieved  reliably. 
Progressive  improvements  in  the  hardware  and  instrumentation 
of  the  turbine  test  rig  following  each  series  of  tests,  and 
finally  the  simplification  of  the  data  reduction  program, 
resulted  in  the  successful  separation  of  losses  for  particular 
configurations  reported  in  Reference  5.   However,  anomalies 
remained  to  be  explained ,  and  the  values  obtained  for  the 
stator  and  rotor  loss  coefficients  were  scattered. 

The  questions  raised  by  the  preceding  test  program  mainly 
involved  the  results  of  the  turbine  configuration  designated 
in  Reference  5  as  Turbine  C.   Turbine  C  had  converging- 
diverging  stator  passages  in  an  axial  entry,  single  impulse 
stage  which  was  designed  to  operate  at  a  pressure  ratio  of 
3.71.   Results  obtained  for  this  geometry  discharging  to 
atmosphere  did  not  agree  with  results  obtained  at  reduced 
pressure  levels  when  plotted  in  terms  of  similarity  varia- 
bles ("referred"  quantities) .   Furthermore,  the  values 
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obtained  for  flow  rate  were  of  questionable  accuracy  since 
the  existing  flow  nozzle  was  too  large.   Also  the  correction 
for  a  labyrinth  leak  rate  became  a  more  significant  fraction 
of  the  turbine  flow  rate. 

The  goal  of  the  present  work  was  to  determine  the  blade 
row  and  stage  performances  of  Turbine  C,  while  resolving  the 
anomalies  reported  in  Reference  5.   This  would  allow  the  use 
of  the  test  rig  to  investigate  effects  of  geometrical  changes 
(stator-rotor  separation  and  tip  clearance  for  example)  on 
the  losses.   The  first  step  was  to  design  and  construct  a  new 
flow  nozzle  for  the  range  of  flow  rates  required  by  Turbine 
C.   The  flow  nozzle  was  carefully  calibrated  against  three 
different  nozzles  operating  choked.   The  leakage  rate  for 
the  labyrinth  seal  was  also  measured  and  a  semi-analytic 
representation  of  the  leakage  rate  as  a  function  of  pressure 
ratio  was  obtained  which  should  apply  at  different  tempera- 
tures.  These  developments  are  reported  in  Appendix  A.   The 
existing  data  reduction  program  [Ref.  5]  was  revised  to 
access  a  Hewlett-Packard  Model  9  867B  Mass  Memory  unit.   The 
revised  program  is  described  fully  in  Appendix  B. 

In  the  results  of  the  test  program  given  here  in  Section 
III,  the  effect  of  pressure  level  reported  in  Reference  5 
was  not  observed.   The  scatter  in  the  loss  coefficient  however 
was  not  removed  by  the  improvements  in  the  accuracy  of  the 
measurements.   In  the  discussion  in  Section  IV,  a  method 
for  smoothing  the  loss  coefficient  data  is  demonstrated 
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which  will  then  allow  a  useful  study  of  the  effect  of  axial 
and  tip  clearances  to  be  carried  out. 
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II.   TURBINE  TEST  RIG  INSTALLATION 

A.   DESCRIPTION 

The  test  installation  consists  of  three  major  components: 
an  Allis  Chalmers  twelve  stage  axial  flow  compressor,  an 
exhauster  assembly,  and  the  turbine  test  rig  (TTR)  itself. 

The  compressor  is  the  source  of  driving  air  for  the  TTR 
and  for  the  exhauster  assembly.   Fig.  A-l  shows  the  piping 
arrangement.   Turbine  air  passes  through  the  first  settling 
tank  into  an  eight-inch  pipe  containing  a  flow  nozzle,  into 
the  second  settling  tank  and  into  the  turbine. 

Fig.  1  shows  the  plenum,  the  floating  stator  assembly, 
the  rotor,  and  the  dynamometer  [Ref.  5],   Pressure -ratios  of 
6:1  can  be  achieved  when  the  system  is  hooded.   The  hood 
was  needed  to  achieve  near  design  pressure  ratios  in  the 
tests  reported  here.   Fig.  2  shows  the  turbine  blading 
arrangement  of  the  stator  and  rotor.   Ref.  8  contains  detailed 
description  of  the  test  rig  hardware. 

The  floating  stator  assembly  depicted  in  Fig.  3  permits 
measurements  of  the  stator  torque  while  axial  and  rotational 
movements  are  detected  by  calibrated  force  transducers  that 
are  heat  insensitive.   These  measurements  allow  the  deter- 
mination of  the  axial  and  tangential  velocity  components 
at  the  stator  exit  [Appendix  A,  Ref.  5]. 

In  this  report  one  configuration  designated  Turbine  C 
was  tested,  the  geometry  of  which  is  shown  in  Fig.  4;  Table 
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1  describes  the  geometry  quantitatively.   The  stator  blade 
profile  is  shown  in  Fig.  5.   The  blades  generate  a  converging- 
diverging  nozzle  shape.   Pressure  measurements  are  taken  at 
the  locations  shown  in  Fig.  5.   The  pressures  necessary  to 
the  analysis  of  stator  forces  (Appendix  B)  are  taken  at  the 
locations  shown  in  Fig.  6 

B.  TEST  MEASUREMENTS  AND  ACCURACY 

MASS  FLOW  RATES :  Appendix  A  gives  a  detailed  analysis 
of  the  turbine  flow  rate  and  the  labyrinth  seal  leak  rate, 
and  the  results  of  calibration  measurements.  The  results 
for  the  discharge  coefficient  are  given  in  Fig.  A-6.  Fig. 
A-10  shows  the  Kinetic  Energy  Factor  (Appendix  A)  obtained 
in  labyrinth  seal  leak  rate  calibration  and  its  respective 
polynomial.  The  equations  used  in  calculating  the  flow  rate 
from  measurements  are  summarized  in  Table  A-l  and  Table  A-2. 

FORCES,  TORQUES,  TEMPERATURES  AND  PRESSURE:   References 
8  and  9  give  calibration  procedures  for  the  TTR.   Identical 
procedures  were  employed  here.   Table  II  of  Ref.  5  gives  the 
expected  accuracies  of  these  measurements. 

C.  TESTING  AND  DATA  REDUCTION 

The  TTR  data  collection  system  is  described  in  Ref.  11. 
Appendix  D  in  Ref.  5  gives  a  detailed  explanation  of  the 
turbine  test  procedures.   The  data  reduction  program  developed 
in  Ref.  5  was  revised  to  use  the  mass  memory  system.   The 
revised  program  is  described  in  Appendix  B. 
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Seven  tests  were  conducted  of  Turbine  C  at  pressure 
ratios  from  1.75  to  3.25.   The  parameters  and  conditions 
of  the  tests  are  summarized  in  Table  II. 
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III.       RESULTS 

A.  INTRODUCTION 

The  reduced  data  from  seven  tests  of  Turbine  C  are 
listed  in  Table  V  to  XI.   Table  IV  gives  a  listing  of  the 
symbols  and  meanings  of  the  Table  headings  used  in  Table  V 
to  XI.   Plots  of  the  reduced  data  are  presented  in  Figures 
7  to  16.   Table  III  gives  an  explanation  of  symbols  used  in 
these  figures. 

Tests  1  and  2  were  conducted  with  the  hood  off.   Tests 
3,  4,  5,  6  and  7  were  made  with  the  hood  installed.   Tests 
5  and  6  were  repeats  of  tests  3  and  4  respectively.   The 
results  are  discussed  in  the  following  paragraphs.  ' 

B.  TORQUE 

Fig.  7a  and  7b  give  the  results  of  the  referred  torque 
versus  referred  speeds  at  a  constant  pressure  ratio.   A 
nearly  linear  dependency  exists  for  all  seven  runs  with 
consistency  in  performance  for  the  1.75  and  2.25  hooded  and 
unhooded  pressure  ratios. 

The  results  from  the  unhooded  tests  (Fig.  7a)  agree  with 
the  results  of  the  hooded  tests  (Fig.  7b) ,  at  the  same 
pressure  ratio,  to  within  1  percent. 

C.  EFFICIENCY  (TOTAL-TO-STATIC) 

Figures  8a  and  8b  give  the  results  of  the  total-to-static 
efficiency  versus  referred  R.P.M.   Figures  9a  and  9b  show 
plots  of  efficiency  versus  isentropic  head.   The  maximum 
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efficiency  occurred  at  similar  values  of  the  isentropic 
head  for  all  pressure  ratios.   Ref.  17  relates  efficiency 
to  isentropic  head.   The  peak  efficiency  was  achieved  at  an 
isentropic  head  coefficient  between  4.0  and  4.5  for  both  the 
hooded  and  unhooded  tests. 

The  efficiency  is  seen  to  increase  as  the  pressure 
ratio  is  increased,  approaching  80  percent  at  the  pressure 
ratio  of  3.25.   The  efficiency  would  be  expected  to  increase 
to  the  design  pressure  ratio  of  3.71,  at  which  the  stator 
nozzles  should  be  just  correctly  expanded. 

D.   THEORETICAL  DEGREE  OF  REACTION 

Figures  10a  and  10b  give  the  results  of  theoretical 
degree  of  reaction  versus  isentropic  head  coefficient  at 
different  pressure  ratios.   Ref.  5  and  Ref.  17  give  the 
equation  for  the  theoretical  degree  of  reaction.   For  all 
runs  the  theoretical  degree  of  reaction  was  negative  at 
lower  speeds.   For  a  negative  value,  P2/P.   must  be  greater 
than  P./P   .   The  values  become  positive  only  at  high  r.p.m. 
in  both  the  hooded  and  unhooded  configurations,  passing 
through  the  design  value  of  zero  at  a  value  of  the  isentropic 
head  coefficient  corresponding  to  the  maximum  in  efficiency. 
Thus,  negative  values  of  the  theoretical  degree  of  reaction 
occur  at  speeds  lower  than  the  design  speed  at  a  given 
pressure  ratio.   Typically,  there  is  seen  to  be  also  a  loss 
of  three  to  five  percent  in  efficiency  from  the  maximum 
value  as  the  isentropic  head  coefficient  is  increased  100 
percent  above  the  design  value. 
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Results  from  hooded  and  unhooded  tests  are  seen  to 
agree  with  the  exception  of  Run  3  in  Figure  10b  for  which 
the  values  are  low.   Lack  of  agreement  for  the  data  of  this 
test  was  found  in  all  but  the  overall  stage  performance 
measurements . 

E.  EFFECTIVE  DEGREE  OF  REACTION 

Figures  11a  and  lib  give  the  results  for  the  effective 
degree  of  reaction  versus  isentropic  head  coefficient.   At 
all  pressure  ratios,  hooded  and  unhooded,  the  effective 
degree  of  reaction  was  negative  at  lower  speeds.   Ref.  17 
states  that  the  effective  degree  of  reaction  is  a  means  for 
judging  whether  the  flow  is  being  accelerated  or  decelerated 
in  the  rotating  cascade.   The  negative  degree  of  reaction 
in  Fig.  11a  and  lib  show  that  the  flow  is  decelerated  at 
lower  r.p.m.  and  becomes  accelerated  as  r.p.m.  is  increased 
(Isentropic  head  is  decreased) .   In  no  case  should  the  design 
value  of  the  effective  degree  of  reaction  be  less  than  zero. 
As  stated  above,  operation  at  speeds  less  than  design  speed 
results  in  a  loss  in  stage  efficiency,  partially  due  to  the 
undesirable  deceleration  in  the  rotor. 

The  results  for  Run  4  (at  a  pressure  ratio  of  3.2  5)  in 
Figure  lib  reflect  the  marked  improvement  in  efficiency  at 
pressure  ratios  approaching  design. 

F.  STATOR  LOSS  COEFFICIENT 

The  results  for  the  stator  loss  coefficient  are  shown 
in  Fig.  12a  and  Fig.  12b.   While  there  is  some  scatter  in 
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the  data,  the  results  from  hooded  and  unhooded  tests  are 
seen  to  be  in  reasonable  agreement.   There  is  a  slight 
increase  in  the  stator  loss  as  speed  is  increased  at  any 
pressure  ratio.   The  stator  loss  for  the  highest  pressure 
ratio  (3.25)  is  significantly  lower  (13-17  percent)  than  for 
the  lower  pressure  ratios  (18-25  percent).   This  reflects 
the  poor  performance  of  a  converging-diverging  nozzle  at 
pressure  ratios  less  than  design. 

G.   ROTOR  LOSS  COEFFICIENT 

The  results  for  rotor  loss  coefficient  are  shown  in 
Figure  13a  and  Figure  13b.   For  completeness,  the  results  of 
Run  3  have  been  included  where  the  rotor  loss  coefficients 
were  evaluated  to  be  positive.   Negative  loss  coefficients 
were  obtained  at  higher  speeds  in  Run  3,  indicating  an 
error  in  at  least  one  element  of  the  primary  data. 

The  rotor  loss  coefficient  is  the  performance  parameter 
most  sensitive  to  small  errors  (and  possibly  unsteadiness) 
in  the  measurements.   The  sensitivity  arises  because  it 
involves  most  of  the  parameters  (see  Ref.  5)  previously 
calculated  in  the  form  of  a  ratio  of  small  differences.   In 
view  of  this,  if  the  results  of  Run  3  in  Fig.  13b  are  excluded 
from  consideration,  the  rotor  loss  coefficient  for  the  tests 
are  reasonably  consistent.   The  results  for  hooded  and 
unhooded  tests  agree  to  within  20  percent  at  corresponding 
pressure  ratios.   There  is  also  a  definite  decrease  in  the 
rotor  loss  coefficient  as  pressure  ratio  is  increased  toward 
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the  design  value,  and  a  consistent  decrease  as  the  speed 
is  increased  at  a  fixed  pressure  ratio.   The  simultaneous 
reduction  observed  in  the  rotor  and  stator  loss  coefficients 
at  higher  pressure  ratios  is  consistent  with  the  increase 
measured  in  the  stage  efficiency  (Fig.  8  and  Fig.  9) . 

H.   FLOW  ANGLES 

An  example  of  the  variation  in  the  flow  angles  at  the 
stator  and  rotor  exits  as  a  function  of  the  isentropic 
head  coefficient  is  shown  in  Fig.  14.   As  was  found  in  all 
cases  tested  in  Ref.  5,  the  variations  in  stator  exit  flow 
angle  and  rotor  exit  relative  flow  angle  are  very  small. 
Also,  the  total  variation  of  the  relative  flow  angle  into 
the  rotor  is  only  about  6  degrees  for  a  50  percent  change 
in  r.p.m.   The  results  in  Fig.  14  are  qualitatively  repre- 
sentative of  the  results  for  each  of  the  seven  tests.   The 
exact  variations  in  the  measured  angles  can  be  seen  in 
Tables  V  to  XI . 

The  stator  exit  flow  angle  is  shown  in  Fig.  15a  and 
Fig.  15b  as  a  function  of  isentropic  head  coefficient  for 
the  seven  tests.   It  can  be  seen  that  the  measurements  are 
in  the  range  of  75  to  77  degrees,  with  the  lower  value 
occurring  at  the  highest  pressure  ratio  (3.25).   It  can  be 
seen  in  Fig.  5  that  the  suction  surface  of  the  stator  blade 
is  designed  to  exit  at  75  degrees,  and  the  pressure  surface 
at  a  smaller  angle.   Thus  the  measured  stator  exit  angle 
appears  to  be  slightly  larger  than  would  be  expected;  however, 
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the  measured  values  when  the  nozzle  passages  are  nearly 
correctly  expanded  approach  the  geometrical  angle. 

The  large  values  of  stator  exit  flow  angle  at  lower 
pressure  ratios,  and  the  higher  values  of  stator  loss 
coefficient  are  probably  the  results  of  the  presence  of  flow 
separation.   The  effect  of  flow  separation  would  be  reflected 
in  two  ways.   Firstly,  the  flow  angle  of  the  separated  flow 
could  depart  significantly  from  the  trailing-edge  angle. 
Secondly,  the  existence  of  a  region  of  separated  flow  could 
seriously  affect  the  averaging  process  which  is  implicit 
in  the  data  reduction. 

I.   FLOW  RATE 

The  referred  flow  rate  is  shown  as  a  function  of  the 
isentropic  head  coefficient  for  the  seven  tests  in  Fig. 
16a  and  Fig.  16b.   The  referred  flow  rate  was  constant  at 
1.02  lbs/sec  to  within  the  accuracy  of  the  measurements 
for  all  speeds  and  all  pressure  ratios  when  operating 
hooded  (Fig.  16b).   During  unhooded  operation  (Fig.  16a), 
the  results  gave  a  constant  value  of  1.01  lbs/sec. 


22 


IV.   DISCUSSION 

A.  RESOLUTION  OF  PREVIOUS  ANOMALIES 

In  the  results  presented  above,  the  difference  in  the 
turbine  performance  operating  hooded  compared  with  unhooded, 
is  seen  to  be  small.   This  is  in  contrast  to  the  results  pre- 
sented in  Reference  5.   Three  factors  contributed  to  the 
resolution  of  this  anomaly;  the  accurate  measurement  of  the 
turbine  flow  rate,  installation  of  temperature  insensitive 
load  cells  within  the  hood  to  measure  stator  loads,  and 
corrections  that  were  made  to  the  data  reduction  program. 
There  remains  to  be  explained  a  difference  of  1%  in  the 
flow  rate  between  hooded  and  unhooded  operation.   More  data 
are  needed  from  unhooded  tests  to  confirm  that  this  differ- 
ence is  repeatable.   Because  of  the  sensitivity  of  the  perfor- 
mance parameters  to  the  flow  rate,  the  residual  disagreement 
between  hooded  and  unhooded  results  might  be  resolved  if 
the  difference  in  flow  rates  were  explained. 

B.  FLOW  RATE 

The  use  of  a  choked  nozzle  to  calibrate  a  large  flow 
nozzle  was  successful  (Appendix  A) .   With  the  blockage 
factor  correction  applied,  the  results  for  different  sizes 
of  the  choked  nozzle  agreed  in  the  range  for  which  the  flow 
rates  overlapped. 

The  labyrinth  leak  rate  calibration  measurements  were 
well  represented  by  the  analysis  given  in  Reference  7.   The 
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value  obtained  for  the  Kinetic  Energy  factor  was  higher  than 
values  given  in  Reference  5;  however,  the  geometry  of  the 
turbine  testing  labyrinth  was  not  known  exactly.   There  was 
some  eccentricity  in  the  alignment  of  the  stator  section 
within  the  gland,  so  that  the  clearance  was  nonuniform  to 
an  unknown  degree . 

The  flow  rate  through  the  turbine,  when  expressed  as 
a  referred  quantity,  was  measured  to  be  a  constant  over  the 
range  of  pressure  ratios  tested.   (However,  there  was  scatter 
in  the  measurements  that  resulted  from  unsteadiness  in  the 
measurements  of  nozzle  pressure  drop.)   Since  the  stator 
was  choked,  the  constant  value  was  expected.   The  value 
of  the  constant  was  consistent  with  the  flow  rate  computed 
for  the  stator  using  a  calculated  value  for  the  stator 
throat  area.   A  value  of  the  blockage  factor  within  1%  of 
unity  was  required  to  be  applied  to  the  geometrical  stator 
area  to  obtain  total  agreement.   Since  the  exact  area  of  the 
throat  where  choking  occurs  in  the  stator  is  not  known  to 
this  accuracy,  it  is  concluded  that  flow  nozzle  and  stator 
throat  measurements  of  the  flow  rate  are  in  agreement. 

C.   LOSS  COEFFICIENTS 

The  scatter  observed  in  the  results  for  the  loss  coeffi- 
cients makes  a  comparison  between  predicted  and  measured 
values  inconclusive.   There  are  trends  in  the  data  which  are 
discernible,  the  rotor  losses  are  consistently  higher  than 
the  stator  losses,  for  example.   However,  the  scatter  does 
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not  allow  the  effect  of  parameter  changes  on  the  separate 
loss  coefficients  to  be  determined.   The  scatter  is  the 
result  of  the  extreme  sensitivity  of  the  loss  coefficient  to 
the  separate  measurements  on  which  they  depend.   For 
example,  the  stator  loss  coefficient  (z;  )  is  given  by 

s.  =  1 ~ET  (1) 

where  X,  is  the  non-dimensional  velocity  and  p,  is  the 
calculated  pressure  ratio  (P-./P.  )  at  station  1.   Differ- 
entiation of  Eq.  (1)  gives,  if  X-,  is  constant, 


III 

d^S         ^S   Y-l      Pl  Y  dPl 

2.  =  -r( E.)  (lZ±.)  ( i ±  (7) 

Cs      U  Cs  M  Y  M       Y-D  3  PX            U) 

1  -  P.  y 


Taking  the  sonic  condition  of  P,  =  0.532,  with  y   =   1 . 4  and 
r   =  0.2,  the  quantity  in  brackets  has  the  value  5.78.   Thus 
a  1%  error  in  calculating  P,  produces  -6%  error  in  the  stator 
loss  coefficient.   Since  P,  is  derived  (as  is  X,)  from 
measurements  scanned  over  a  period  of  more  than  a  minute, 
with  some  variation  in  operating  point  (as  well  as  an 
observable  fluctuation  in  the  flow  nozzle  pressure  drop) 
occurring  during  the  data  scan,  the  scatter  in  the  loss 
coefficient  data  is  understandable. 
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Before  the  loss  coefficient  can  be  used  as  a  measure 

of  performance,  a  method  of  smoothing  the  data  while 

retaining  accuracy  is  needed.   Two  attempts  were  made 

here.   First,  an  attempt  was  made  to  smooth  the  calculated 

values  of  the  interstage  pressure  (p, ) ,  by  comparing  p, 

to  the  pressures  measured  at  the  hub  and  tip.   Figure  17 

shows  the  variation  p,  and  the  measured  pressure  with 

isentropic  head  coefficient  for  one  particular  test. 

Clearly,  the  calculated  average  pressure  at  station  1  must 

be  smoothly  behaved  if  the  hub  and  the  tip  pressures  are 

smoothly  behaved.   An  attempt  to  derive  a  smoothing  function 

for  the  interstage  pressure  is  shown  in  Fig.  18,  where  the 

data  from  several  tests  at  the  same  pressure  ratio  are 

plotted  in  a  dimensionless  form.   The  figure  clearly  shows 

the  degree  of  scatter  in  the  values  of  P, .   On  the  basis  of 

these  data  above,  however,  the  use  of  a  single  polynomial 

for  the  quantity  (P,-P,  1_)/(PJ_   -P.  ,)  as  a  function  of 
^      J         1   hub    tip   hub 

isentropic  head  coefficient  cannot  be  justified.   When  a 
polynomial  was  derived  from  the  data  of  a  single  test,  and 
the  loss  coefficients  for  the  stator  and  rotor  were  recomputed, 
the  results  shown  in  Fig.  (19a)  and  (19b)  were  obtained. 

The  second  attempt  to  smooth  the  loss  coefficient 
data  involved  effectively  the  elimination  of  scatter  in  the 
input  data  from  which  the  interstage  pressure  and  velocity 
were  calculated.   The  referred  flow  rate  was  determined  to 
be  constant  when  the  stator  was  choked.   By  setting  the 
referred  flow  rate  equal  to  the  constant  value  (1.02  or  1.01, 
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hooded  and  unhooded,  respectively)  and  continuing  the  data 
reduction  as  before,  the  results  shown  in  Fig.  2  0a  and  20b 
for  the  loss  coefficient  were  obtained.  A  considerable 
reduction  in  scatter  is  evident,  with  the  major  benefit  in 
the  rotor  loss  coefficient.  Further  data  must  be  examined 
to  fully  appraise  the  benefit  of  this  smoothing  technique. 

A  third  technique,  so  far  not  attempted,  is  to  calcu- 
late P,  in  the  data  reduction  before  calculating  the  axial 
velocity  component  at  station  1.   The  values  of  P,  could 
be  smoothed  using  a  polynomial  curve  fit  and  the  smoothed 
values  of  P,  used  to  calculate  the  velocity.   An  application 
of  all  three  techniques  should  produce  loss  coefficient 
distributions  which  are  smooth  and  yet  accurate. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  are  made  for  turbine  configura- 
tion C. 

1.  The  loss  coefficients  are  extremely  sensitive  to 
input  data,  especially  the  rotor  loss  coefficient.   Repeata- 
bility in  calculating  the  loss  coefficient  directly  from 
measurements  was  poor. 

2.  Methods  were  demonstrated  to  smooth  the  loss  coef- 
ficient data  without  sacrificing,  and  possibly  improving, 
the  accuracy. 

3.  For  all  pressure  ratios  there  was  a  velocity 
deceleration  through  the  rotor  and  a  pressure  rise  'across 
the  rotor  at  lower  r.p.m.,  indicating  perhaps  a  shock  wave 
at  or  in  the  rotor. 

The  following  recommendations  are  made: 

1.  Apply  the  three  data  smoothing  techniques  described 
above  to  all  data  from  Turbine  C  in  sequence:   input  constant 
referred  flow  rate,  calculate  interstage  pressure  and  smooth 
using  measured  hub  and  tip  pressures,  then  calculate 
interstage  velocity. 

2.  Determine  a  more  accurate  input  for  K, ,  the  ratio 

of  the  momentum  average  velocity  to  the  mass-average  velocity 
and  examine  the  effect  on  the  calculated  stator  exit  flow 
angle. 
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3.   Program  the  quadratic  curves  for  the  theoretical 
loss  coefficient  in  references  16  and  18  in  order  to  reduce 
the  input  error. 

Finally,  it  is  concluded  that  following  the  recommenda- 
tions listed  here,  the  effect  of  parameter  variations  on 
the  blade  row  losses  could  be  examined  satisfactorily  using 
the  turbine  test  rig. 
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TABLE  I 
TURBINE  GEOMETRY 


TURBINE 

STATOR 

ROTOR 

C 

1 

1 

DESCRIPTION 

A*   2 
(in? 

NUMBER 
OF 
BLADES 

\ix) 

(in.) 

STATOR  1 

CONVERGING 
DIVERGING 
NOZZLES 

2.9058 

31 

4.184 

0.5775 

ROTOR  1 

CIRCULAR 
ARC 

7.119 

60 

Rm1=4.193 
Rm2=4.250 
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h2=0.8475 

30 


w 

S3 

H 

PQ 
04 


<  en 

ft  w  z 

H  hi  H 

H   e_5  ^         O 


O 

i— i 

o 


o 

I— I 

o 
o 


o 

r- 1 

o 
o 


o 

I— I 
o 

o 


o 

.—I 

o 
o 


o 

1— I 

o 

o 


o 

I— I 

o 
© 


J  M  /-»  O  O 

<!<![/]  CM  CM 

HW2  CM  CM 

<  o  ^  o  o 


o 

CM 
CM 


o 

CM 
CM 


o 

CM 
CM 


o 

CM 
CM 


o 

CM 
CM 


w 

ptf 

S  CM 

tn  O  ft 

63  13     ° 

ft  <      -u 
ft  ft  ft 


<|- 

CM 


CO 
■U 
CO 

c 


T3 
O 
O 

4J 


CO 
■U 

o 

3 

fl 

O 
O 

<L> 

U 

cu 

5 


m 


CO 

C 

3 
ft 


CM 


m 

CM 


m 


m 

CM 


m 

CM 


CU 

C 

•r-l 

.D 
U 

3 

<4-l 

O 

>, 
U 
■U 

cu 

e 
o 
cu 
Ml 

CO 

cu 
> 

•H 


CO 

H 


CM 


CM 


co 


CM 


CM 


>» 

CJ 

a 

cu 

•H 

CM 

o 

ft 

4-1 

o 

14-4 

4J 

w 

ft 

v_^ 

^> 

II 

o 

.— 1 

O 

•H 

•H 

4-1 

o 

4J 

CO 

H 

CO 

4-1 

ft 

C/2 

r-- 

1 

CU 

o 

c/) 

o 

h 

■u 

W 

s 

3 

i 

ft 

M 

CO 

i— i 

3 

g 

CO 

CO 

O 

5 

cu 

4-1 

HH 

H 

J-l 

o 

ft 

s 

ft 

H 

en 

l-J 

O 
PQ 


Cn 

ft 

o 

O 
H 
H 
< 


t-4 
ft 
X 

m 


o 

PQ 


en 


ft 
ft 


en 

H 


ft" 


CM 


CO 


m 


v£> 


H 
O 


ft 


31 


TABLE  IV 


DEFINITION  OF  COLUMN  HEADINGS  IN  TABLES  V  TO  XI 

SYMBOL  DEFINITION 

VI  VELOCITY  AT  STATOR  EXIT  PLANE  (FT/SEC) 

V2  VELOCITY  AT  ROTOR  EXIT  PLANE  (FT/SEC) 

VA1  AXIAL  VELOCITY  AT  STATOR  EXIT  PLANE  (FT/SEC) 

VA2  AXIAL  VELOCITY  AT  ROTOR  EXIT  PLANE  (FT /SEC) 

VUI  TANGENTIAL  VELOCITY  AT  STATOR  EXIT  PLANE  (FT/SEC) 

VU2  TANGENTIAL  VELOCITY  AT  ROTOR  EXIT  PLANE  (FT/SEC) 

MI  MACH  NUMBER  AT  STATOR  EXIT  PLANE 

MAI  AXIAL  MACH  NUMBER  AT  STATOR  EXIT  PLANE 

M2  MACH  NUMBER  AT  ROTOR  EXIT  PLANE 

MA2  AXIAL  MACH  NUMBER  AT  ROTOR  EXIT  PLANE 

Al  FLOW  ANGLE  AT  STATOR  EXIT  PLANE  (DEGREES) 

A2  FLOW  ANGLE  AT  THE  ROTOR  EXIT  PLANE  (DEGREES) 

BI  RELATIVE  FLOW  ANGLE  AT  ROTOR  INLET  PLANE  (DEGREES) 

B2  RELATIVE  FLOW  ANGLE  AT  ROTOR  EXIT  PLANE  (DEGREES) 

ZS  STATOR  LOSS  COEFFICIENT 

ZSTH  THEORETICAL  STATOR  LOSS  COEFFICIENT 

ZR  ROTOR  LOSS  COEFFICIENT 

ZRTH  THEORETICAL  ROTOR  LOSS  COEFFICIENT 

ZR*  ROTOR  CARRY  OVER  LOSS  COEFFICIENT 

ZI  ROTOR  INCIDENCE  LOSS  COEFFICIENT 

P.R.  PRESSURE  RATIO  P   /P 

to   2 

K-IS  ISENTROPIC  HEAD  COEFFICIENT 

NREF  REFERRED  RPM  (R.P.M.) 
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WREF 

ETA 

HPREF 

RTH 

REFF 

PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

P9 

KBLOCK 

STPR 

RPM 


REFERRED  FLOW  RATE  (  LBM/SEC  ) 

TOTAL  TO  STATIC  EFFICIENCY 

REFERRED  HORSEPOWER  (  H.P.  ) 

THEORETICAL  DEGREE  OF  REACTION 

EFFECTIVE  DEGREE  OF  REACTION 

PRESSURE  RATIOS  CORRESPONDING  TO  TAP  LOCATIONS 

IN  FIGURE  5.  PRESSURES  ARE  REFERRED  TO  STATOR 

INLET  TOTAL  PRESSURE 


STATOR  BOCKKAGE  FACTOR 


PRESSURE  RATIO  ACROSS  THE  STATOR  P   /P, 

to   1 


ROTOR  SPEED  (  R.P.M.  ) 
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APPENDIX  A 
FLOW  RATE  DETERMINATION 

A-l   INTRODUCTION 

A  precise  measurement  of  the  flow  rate  is  necessary 
if  the  performance  of  an  individual  blade  row  is  to  be 
determined  accurately  from  measurements  made  with  the  turbine 
test  rig.   Anomalies  that  occurred  in  the  results  from  the 
first  tests  of  the  present  supersonic  turbine  (Turbine  C 
in  Ref.  5)  were  possibly  due  to  inaccuracies  in  the  speci- 
fied flow  rates.   Since  the  flow  rate  for  Turbine  C  was 
smaller  than  for  the  turbines  previously  tested  in  the  rig, 
small  pressure  differentials  were  measured  with  the  existing 
nozzle.   Also,  the  higher  pressure  ratios  at  which  Turbine 
C  was  operated  resulted  in  increased  leakage  through  the 
labyrinth  seals ,  the  geometry  of  which  might  have  change 
slightly  on  reassembly.   Therefore,  before  new  tests  of 
Turbine  C  were  begun,  a  new  flow  nozzle  was  provided  and 
careful  measurements  were  made  of  the  labyrinth  leak  rate. 
The  flow  nozzle  design  and  calibration  are  described  in 
Section  A-2.   The  labyrinth  leak  rate  determination  is  given 
in  Section  A- 3. 

A-2   MEASUREMENTS  OF  THE  TOTAL  FLOW  RATE 

A-2.1   FLOW  NOZZLE  DESIGN 

The  flow  nozzle  is  positioned  in  a  pipe  through  which 
air  is  supplied  to  the  turbine  (Fig.  A-l) .   The  nozzle 
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generates  a  differential  pressure  to  which  the  flow  rate 
is  related. 

For  a  one-dimensional  steady  flow,  the  mass  flow  rate, 
W,  is  given  by 

W  =  p  AV  (A-l) 

where  p  is  the  density,  V  is  the  velocity,  and  A  is  the 
cross  sectional  area. 

Using  the  perfect  gas  equation  of  state 


pt 
p  =  ss-  (A-2) 

t 


and  Bernoulli  equation 


Pt  "  P  =  \   P  V2  (A-3) 


Eq.  (A-l)  for  low  velocities  becomes 


(A-4) 


The  symbols  P  and  T  are  pressure  and  temperature,  respec- 
tively, and  the  subscript  t  denotes  stagnation  conditions. 
This  equation,  with  explicit  corrections  for  compressibility 
and  thermal  expansion,  is  used  to  calculate  the  flow  rate 
through  a  flow  nozzle  from  measurements  of  temperature  and 
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pressures.   Withe  the  addition  of  the  correction  factors  and 
with  conversion  factors  added  to  account  for  the  units  of 
measurements  [Ref .  1] ,  Eq.  (A- 4)  becomes 


0.16384     _  2  _  v      v   /PN  Ah  ...   ,     .    ,_  _, 
W  = D   A^  K^  Y   /— 5—  (lbs/sec.)    (A-5) 

/2.036      N    w   w  x  J   lt 


where 


DN  =  Diameter  of  nozzle  (inches) 


A.,  =  Thermal  expansion  coefficient 
Y.  =  Compressibility  coefficient 


PN  =  Pressure  at  the  nozzle  (psia) 


Ah  =  Water  differential  (in.  H20) 


T.  =  Temperature  (°R) 
0.16384 


=  Conversion  factor 


/2.036 


and  lastly,  KN  is  the  discharge  coefficient. 

The  discharge  coefficient  accounts  for  the  fact  that  the 
flow  is  not  one-dimensional,  and  that  the  pressures  measured 
are  not  P.  and  (P.  -  P) .   Generally,  the  discharge  coefficient 
has  a  value  above  0.95  for  a  standard  ASME  nozzle  at  high 
Reynolds  numbers.   Establishing  the  relationship  between 
Reynolds  number  and  the  discharge  coefficient  is  what  con- 
stitutes the  calibration. 
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The  flow  rates  expected  in  the  supersonic  turbines  were 
1.3  to  3.5  lbm/sec.   A  flow  nozzle  with  a  diameter  of  3.25 
inches  was  chosen  (by  applying  Eq.  (A-5)  with  unity  for  the 
coefficients)  to  give  sutiable  water  column  measurements  of 
approximately  9.8  to  58.1  inches.   The  nozzle  design  followed 
the  criteria  given  in  Ref.  1  for  "Low  Beta  series"  designs. 
The  geometry  of  the  flow  nozzle  is  shown  in  Fig.  A-2  and 
views  are  given  in  Fig.  A-3.   The  locations  of  the  pressure 
taps  are  shown  in  Fig.  A-2.   Four  throat  taps  are  located 
1.5  nozzle  diameters  from  the  nozzle  face  and  are  manifolded 
together.   The  upstream  pressure  is  measured  using  taps  in 
the  upstream  flange. 

A-2. 2   FLOW  NOZZLE  CALIBRATION 

The  discharge  coefficient  is  given  in  terms  of  measure- 
ments, if  the  flow  rate  is  determined  by  independent  means. 
One  of  the  most  accurate  methods  of  determining  the  flow 
rate  of  a  gas  is  by  choking  the  flow  at  a  known  area.   This 
is  because  the  mass  flux  is  nearly  uniform  where  choking 
occurs ,  and  because  the  boundary  layer  can  be  kept  thin  in 
a  rapid  smooth  contraction  from  a  relatively  large  pipe. 
Hence,  in  order  to  determine  the  discharge  coefficient  of  the 
new  flow  nozzle  as  a  function  of  Reynolds  number,  the 
arrangement  shown  in  Fig.  A-4  was  used. 
A-2. 2.1  Method 

Flow  nozzle  calibration  runs  were  made  at  con- 
trolled supply  pressures  from  15  psia  to  45  psia.   Three 
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different  orifice  plates  with  diameters  of  1.6,  2.065,  and 
2.24  inches  were  mounted  in  turn  to  provide  choked  flows 
over  the  desired  range  of  Reynolds  numbers.   The  differen- 
tial pressure  at  the  flow  nozzle  was  varied  in  increments 
of  2  to  4  inches  of  water.   The  range  of  flow  rates  over 
which  the  orifice  plates  were  choked  at  the  supply  pressure 
was  from  1.33  to  3.55  lbm/sec.   The  maximum  pressure  ratio 
obtained  was  2.92. 

A-2.2.2   Analysis 

The  discharge  coefficient  was  given  by  Eq.  (A-5) , 
and  the  Reynolds  number  was  determined  by 


48  W   3 
i%   =  ^—  (A-6) 


7T  D    li 
P 


where  \i   is  the  viscosity  and  3  is  the  ratio  of  the  diameter 
of  the  nozzle  to  the  diameter  of  the  pipe.   The  flow  rate, 
which  was  the  same  for  the  nozzle  as  for  the  choked  plate, 
was  determined  by  the  following  equation 


where 


r  P   A  K 

W„  =  r   t  B (A-7) 

N    /RT7? 


P.  =  pressure  at  choked  nozzle  (psia) 

2 
A  =  area  of  choked  nozzle  (in  ) 
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K_  =  blockage  factor 

R  =  gas  constant  (ft-lbf/lbm  -  °R) 


T  =  Temperature  at  choked  nozzle  (°R) 


The  function  r  is  a  known  function  of  the  specific  heats  and 
is  given  by 


=  p~<&W  (A-8) 


The  blockage  factor  K_  accounts  for  the  boundary- 
layer  and  for  the  axi symmetric  surface;  K„  is  given  by 
[Ref.  3] 


46* 
'B        D 


K„  =  1  -  =£-  (A-9) 


N 


where  the  star  on  D„  signifies  sonic  conditions,  and  6* 

N 

is  the  boundary  layer  displacement  thickness. 

The  boundary  layer  growth  was  estimated  using 
the  method  of  Falkner  and  Garner  given  in  Ref .  6 .   The 
displacement  thickness,  6*,  is  given  by 


n 
HTT 


i--5*_ _  ^n+1  (A-10) 

sl 

Reln+1 
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where  s,  is  the  distance  from  the  stagnation  point  to  the 
choking  station.   H  is  the  form  factor,  and  a  and  n  are 
empirical  constants.   ^  is  a  factor  which  depends  on  the 
pressure  gradient  along  the  surface.   Re,  is  the  Reynolds 
number  based  on  conditions  at  the  throat  and  s, . 

It  is  shown  in  Ref.  6  that  the  value  of  ty   for 
a  similar  contraction  was  0.38.   With  ip  =  0.38,  and  with 
H  =  1.4,  a  =  0.0076  and  n  =  6  for  air,  Eq.  (A-10)  becomes 


m   0.0932174 
6      I   0.1429  Sl  {A   1L) 

K6  , 


The  Reynolds  number  is  given  by 


W  s. 

Re,  =  12  ^  —  (A-12) 

1       Ax  yx 


where  y,  is  the  viscosity  at  the  throat  where  the  tempera- 
ture is  5/6  times  the  stagnation  temperature. 

For  each  operating  condition  for  each  orifice 
plate,  Re,  was  calculated  from  Eq.  (A-12)  using  W  calculated 
from  Eq.  (A-7)  with  K_  =  1.0.   6*  was  then  calculated  using 

a 

Eq.  (A- 11)  and  K_  obtained  from  Eq.  (A- 9) .   W  was  recalcu- 
lated  from  Eq.  (A-7) .   It  was  found  that  the  blockage  factor 
could  be  taken  to  be  constant  with  pressure  for  each  choked 
plate  at  the  following  values: 
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PLATE  # 

THROAT  DIAMETER 

Sl 

(ins. ) 

KB 

1. 

1.6 

3.51 

0.9928 

2. 

2.065 

3.10 

0.9935 

3. 

2.24 

3.22 

0.9901 

The  flow  chart  for  the  analysis  of  the  calibra- 
tion test  results  is  given  in  Fig.  A- 5  and  the  equations 
used  are  summarized  in  Table  A-l. 

A-2.2.3   Results 

The  nozzle  calibration  extended  over  a  range 

5  5 

of  Reynolds  numbers  from  0.081  x  10   to  0.213  x  10  .   The 

results  are  given  in  Fig.  A-6,  where  the  data  for  three 

different  orifices  are  shown  for  the  range  of  pressures 

for  which  the  flow  was  choked.   A  second  order  polynomial 

was  found  to  represent  the  flow  nozzle  coefficient  as  a 

function  of  the  Reynolds  number. 


K.,  =  1.0272  -  0.1598  ( — S)  +  0.3805  ( — 5)  2    (A-13) 
*  10*  106 


where  Re   is  given  by  Eq.  (A-6) ,  was  used  to  reduce  data 
from  tests  of  the  turbine  test  rig. 

A-2.3   APPLICATION 

The  discharge  coefficient  is  a  function  of  Reynolds 
number.   However,  in  application  the  flow  rate  is  to  be 
determined  and  therefore  the  Reynolds  number  is  unknown. 
The  determination  of  the  flow  rate  from  measurements  requires 
the  solution  of  three  equations  which  have  the  general  forms: 
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WN  =  ax  KN  (A-14) 


Kg  =  KN  (ReN)  (A-15) 


ReN  =  a2  WN  (A_16) 


These  equations  represent  Eq.  (A-5) ,  Eq.  (A-13) ,  and  Eq. 
(A- 6) ,  respectively,  where  a,  and  a~  are  known  in  terms  of 
the  geometry  and  measurements  for  any  operating  point. 
Clearly,  by  substituting  for  Re   and  K   using  Eq.  (A-15) , 
and  Eq.  (A-16)  in  Eq.  (A-14) ,  a  single  function  for  the 
unknown  W  is  obtained. 

An  iterative  procedure  is  used  to  solve  the  three 
equations  for  the  unknown  flow  rate.  The  procedure  is 
given  in  Appendix  B. 

A- 3   MEASUREMENT  OF  THE  LABYRINTH  LEAK  RATE 

A- 3 . 1  Labyrinth  Leak  Rate  Calibration  Test 
Since  the  leak  rate  of  the  plenum  labyrinth  cannot  be 
measured  during  turbine  operations,  the  relationship  of  the 
leakage  flow  rate  to  supply  and  hood  pressures  and  supply 
temperature  was  determined  in  a  separate  calibration  test. 
For  this  test,  the  main  flow  through  the  turbine  was  blocked 
by  a  plate  and  the  labyrinth  was  allowed  to  leak  into  the 
hood  as  in  normal  operation.   The  flow  was  passed  through  a 
two- inch  pipe  (Fig.  A-4b  &  Fig.  A-l) ,  containing  an  ASME 
standard  sharp-edged  orifice.   The  usual  supply  pipe  was 
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blocked.  The  flow  was  discharged  through  the  labyrinth 
seal  into  the  hood  in  which  the  pressure  was  controlled 
by  the  ejector. 

In  determining  the  leakage  rate  using  the  ASME  orifice, 
vena  contracta  taps  were  used.   The  differential  pressure 
was  measured  in  inches  of  water  on  a  manometer  which  had 
0.1-inch  graduations.   The  upstream  pressures  in  the  lab- 
yrinth and  in  the  nozzle  were  measured  in  inches  of  mercury 
on  U-tube  manometers  with  atmospheric  pressure  as  a  refer- 
ence.  The  supply  pressure  was  varied  from  14  to  83  psia, 
while  the  differential  pressure  across  the  vena  contracta 
taps  varied  from  2.75  to  37.3  inches  of  water.   The  maximum 
pressure  ratio  across  the  labyrinth  was  5.6. 

A- 3 . 2   Calculation  of  the  Leakage  Flow  Rate 
Since  the  flow  orifice  was  in  a  standard  installation 
according  to  the  criteria  given  in  Ref.  1,  Table  5  of  Ref.  1 
was  used  to  obtain  the  discharge  coefficient  as  a  function 
of  Reynolds  number.   In  order  to  program  the  data  reduction, 
the  data  of  Table  5,  for  the  appropriate  values  of  Beta, 
were  first  approximated  in  the  range  of  Reynolds  numbers 
from  8,000  to  1,000,000  by  a  polynomial  as  shown  in  Fig. 
A-7.   The  discharge  coefficient  at  each  operating  point  was 
then  given  by 
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K__  =  0.603698259  +  0.004688755  (—-)  (A-17) 
N                              ReN 

5  5 

-  1.250025  x  10~3  (^—)2   +  2.0320771882  x  10~4  {——)  3 

ReN  ReN 

5  5 

-5   10   4  -7   10   5 

-  1.56980769  x  10  D  (£*—)       +   4.5277098  x  10    (~~) 

ReN  ReN 


The  flow  rate  was  calculated  at  each  point  using  Eq. 
(A-5) ,  using  the  expressions  for  the  compressibility  and 
thermal  expansion  coefficients  given  in  Table  A-2,  and  the 
discharge  coefficient  given  by  Eq.  (A-17) . 
A- 3 . 3  Analysis  of  the  Results 
A- 3 . 3 . 1   Simple  Method 

A  non-dimensional  flow  rate,  or  flow  function, 
for  the  flow  through  the  labyrinth  can  be  defined  as 


WT  /(R/g)  TT 


During  the  calibration  tests,  the  stagnation 
pressure  (P  )  was  varied  in  steps  for  constant  values  of  the 

Li 

hood  pressure  (P,) .   The  stagnation  temperature  was  approxi- 
mately constant.   The  flow  function  was  calculated  for  each 
operating  point  using  Eq.  (A-18)  and  the  values  plotted  as 
a  function  of  the  pressure  ratio  (P./P.) .   The  results  are 
given  in  Fig.  A- 8. 

A  second  order  polynomial  was  found  to  fit  the 
data  to  an  accuracy  of  ±  4%  as  shown  in  Fig.  A- 8: 
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$  =  0.2721  +  0.0365  (Pj/PL)  -  0.2357  (Pj/PL) 2      (A-19) 


While  Eq.  (A-19)  represents  the  calibration  tests 
results  quite  well,  it  is  an  empirical  result  for  a  limited 
variation  in  test  parameters.   Operation  at  different 
temperatures,  for  example,  could  give  a  different  behavior. 
A-3.3.2   Use  of  the  Kinetic  Energy  Factor 
It  is  shown  in  Ref.  7  that,  under  reasonable 
assumptions,  the  flow  through  a  labyrinth  with  many  teeth 
(N)  can  be  represented  by 

$  =  3  Ke  (A-20) 


where 


1-(P,/P_)2 
*  =  /  N-Ln(Ph/PL)  (A~21) 


and  Ke  is  the  "Kinetic  Energy  Factor".   Here,  the  dependence 
of  the  leak  rate  on  the  upstream  and  downstream  pressures 
(3)  has  been  obtained  analytically.   The  Kinetic  Energy 
Factor  is  a  measure  of  how  much  kinetic  energy  is  recovered 
as  stagnation  pressure  as  the  flow  passes  successive  teeth. 
It  is  expected  to  be  a  function  of  Reynolds  number,  and 
hence  would  include  the  correct  behavior  with  stagnation 
temperature. 
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The  geometry  of  the  labyrinth  is  given  in  Fig. 
A- 9.   Using  the  notation  in  Fig.  A- 9,  the  Reynolds  number 
for  the  leakage  flow  is  given  by 


12  W  C 

ReT  =  — (A-22) 

L    *L   yL 


where 


L  =  2  i  D  C  (A-23) 

and  y   is  the  viscosity  based  on  the  labyrinth  supply 
temperature,  T  . 

Li 

For  the  flow  rates  obtained  in  the  calibration 
tests,  $  was  obtained  from  Eq.  (A-18) .   $  was  calculated 
using  Eq.  (A-21)  and  Ke  determined  using  Eq.  (A-20) .   The 
Reynolds  number  was  obtained  using  Eq.  (A-22)  with  Eq.  (A-23) 
The  Kinetic  Energy  Factor  was  plotted  as  a  function  of  the 
Reynolds  number  and  the  result  is  shown  in  Fig.  A- 10.   A 
fourth  order  polynomial  was  found  to  fit  the  data  to  an 
accuracy  of  ±  2%: 


Ke  =  0.601373889  +  4.36695141  x  10~4  (-—) 

10 


-  2.4530860  x  10  '    (-^|)   +  6.5672233  x  10   x  (-^|) 

10  10 


-  6.086573  x  10  15  (-^1)  4  (A-24) 

10 
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It  is  noted  that  the  use  of  Eq.  (A-24)  and  Eq. 
(A-21)  gives  a  factor  of  two  improvement  in  accuracy  over 
the  use  of  the  purely  empirical  result  given  as  Eq.  (A-19) . 
Also,  the  dependence  of  the  Kinetic  Energy  Factor  on  the 
supply  temperature  is  at  least  partially  included  in  the 
representation  as  a  function  of  Reynolds  number. 

The  flow  chart  for  the  calibration  data  reduction 
is  given  in  Fig.  A-ll  and  a  summary  of  the  equations  is 
given  in  Table  A-2. 

A-3.2   APPLICATION 

In  order  to  use  Eq.  (A-24)  to  calculate  the  leak  rate 
at  a  given  operating  point,  an  iterative  procedure  is 
necessary  since  the  Reynolds  number  is  again  a  function  of 
the  flow  rate.   The  procedure  is  similar  to  that  followed 
in  calculating  the  total  flow  rate  (Section  A-2. 3)  and  is 
described  in  Appendix  B. 
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SUMMARY  OF 

FORMULAS  FOE 

TABLE  A-l 
.  NOZZLE  CALIBRATION  DATA  REDUCTION  PROGRAM 

FLOW 
RATE 

% 

= 

0.16384 

/  2.036   ' 

THEMAL 

EXPANSION 

COEFFICIENT 

\ 

= 

1+0.00193 

1"    *N  -  528     1 

100 

COMPRESSIBILITY 
COEFFICIENT 

Y 

1 

= 
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Ah 

P 

N 

BETA 

e 

= 

D 

N 
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N                P 

7 
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R 
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48  W  @ 
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FUNCTION: 

r 
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H 

2 

y+r 

BLOCKAGE 
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K 
B 
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1 

* 
DN 
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LAYER 
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v 
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0.21364 
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1 

REYNOLDS 
NUMBER:* 

Rel 

= 

W  s1 
12  

Plow 

RATE 

CHOKED 

NOZZLE 

W 

= 

pt  a  k  B  r 

1 1  1 

J       RTt/g 
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SUIMARY  OF   ] 

FORMULAS 

TABLE  A- 2 
FOR  LABYRINTH  CALIBRATION  REDUCTION  PROGRAM 

FLOW 
FUNCTION: 

£ 

.        /    RTt        " 

=    z-030  wl  J  '  g "■  " 

PL  AL 

FLOW 

RATE: 

»L 

0.16384 

.       „  2'         „     „       /    pt4h" 

/  2.036                           *       """^J         Tt 

BETA: 

P 

/  i  -  <3i/ 

-/               PL                                    N     =      10 

i    N-  Lff    (Ph/   PL) 

REYNOLDS 
NUMBERS: 

Rel 

12WLC 

KM 

AREA 
LAB: 

*t 

=    2TTd  p 

KINETIC 

ENERGY 

FACTOR: 

Ke 

.     I 

P 

THERMAL 

EXPANSION 

COEFFICIENT: 

A 
N  • 

=     1  +  0.00193 

"    %  -5281 

100        J 

COMPRESSIBILITY 
COEFFICIENT: 

Yi 

-     1-0.05246      (   0.41  +  0.350*  )      

PN 

BETA: 

fi 

DN 
=     * »        D     =  3.25   ins.,   D„  =  7.975   ins. 

D                                                         p 

P 

VISCOSITY: 

M 

-5 
1.153  x  10 

198.72 

x  0.06333/  TN* 

/TN+1 
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LABYRINTH  SEAL 
TRANSONIC  TURBINE  TEST  RIG 
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FIG.  A-9  PLENUM  LABYRINTH 
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APPENDIX  B 
TURBINE  TEST  RIG  (TTR)  DATA  REDUCTION  AND  PROCESSING 

B-l   INTRODUCTION: 

This  appendix  documents  the  current  data  reduction 
program  which  was  revised  and  updated  from  the  program  used 
in  Ref.  5.   The  present  program  uses  the  "Mass  Memory"  disc 
storage  unit  attached  to  the  Hewlett-Packard  Model  9  830A 
Calculator.   The  equations  developed  in  Ref.  5  apply  to 
the  present  program. 

A  description  of  the  programs  is  given  in  Section 
B-2  and  the  variables  are  listed  and  defined.   Procedures 
for  running  the  program  are  given  in  Section  B-3. 

B-2   DESCRIPTION  OF  PROGRAMS: 

The  reduction  program  is  divided  into  eight  separate 
sub-programs.   Fig.  B-l  shows  the  contents  of  each  sub- 
program.  Table  B-l  gives  the  channel  and  port  numbers  for 
each  of  the  test  measurements  and  the  respective  matrix 
elements  assigned  in  the  data  reduction  program.   Table  B-2 
defines  the  variables  and  equations  (Ref.  5  &  16)  used  in 
the  program.   Table  B-3  is  a  record  or  score  sheet  of  the 
variables  used  in  the  program. 

B-3   PROCEDURE  FOR  DATA  REDUCTION  PROGRAM  USING  MASS  MEMORY: 

1 .   CHECK : 

a.)   Mass  Memory  to  UNLOAD  Position 

b.)   ROM  installed 

c.)   NO  cassettes  in  transport 
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2.  Turn  on  HP  9830A  Calculator 

3.  Turn  on  HP9866A  Printer 

4.  Turn  on  HP  9868A  I/O  Expander 

5.  Turn  on  HP9863A  Tape  Reader 

6.  Turn  on  HP  9  867B  Mass  Memory 

7.  Turn  on  HP  11305A  Controller 

8.  Load  platter  number  TTR  when  the  DOOR  UNLOCKED  light  is 
lit. 

9.  Switch  to  LOAD  position;  wait  30  seconds  until  DRIVE 
READY  light  is  lit. 

10.  Key  in:   GET  "TTR  1:"  Wait  for  completion 

11.  Press  "RUN  EXECUTE." 

12.  The  display  will  read  "ENTER  NEXT  RECORD  #  ON  DATA  FILE." 
Input  your  next  record  number  for  raw  data  storage  and 
press  "EXECUTE." 

13.  The  display  will  read  "TAPE:   1ST  HOLE7-ON  START? :-CONT. " 
By  pressing  "CONTINUE,  EXECUTE"  the  paper  tape  will  be 
read. 

14.  After  the  paper  tape  is  read,  the  display  will  read 
"CORRECTIONS  TO  DATA?  YES=1,  NO=0;"   This  option  permits 
you  to  correct  any  erroneous  corrections.   PRESS  "1 
EXECUTE"  if  you  have  corrections.   Press  "0  EXECUTE" 

if  no  corrections;  go  to  step  19. 

15.  With  Corrections  to  be  made,  the  display  will  read 
"PRESS  PRT  ALL  KEY  FOR  RECORD."   This  enables  you  to 
have  a  record  of  your  data  points. 
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16.  The  display  will  read  "ENTER  CORRECTION  AS  MATRIX 
ELEMENTS."   When  a  scanivalve  is  detected  to  have  given 
an  erroneous  value,  the  manometer  board  can  be  read 
and  correct  value  recorded.   See  Table  B-2  for  the 
appropriate  matrix  identification. 

17.  The  display  will  read  "ENTER  CORRECT  VALUE  =  ?;" 
Input  the  correct  value  and  press  "EXECUTE." 

18.  The  display  will  read  "ANY  MORE  CORRECTIONS?  YES=1, 
NO=0."   Press  "1  EXECUTE"  if  more  corrections;  go  to 
step  17. 

19.  1.5  inches  of  paper  will  emerge  from  the  printer  and 
the  raw  data  will  be  printed.   Afterwards  the  display 
will  read  "STORE  DATE?  ENTER  YES=1,  NO=0."   If.  raw 
data  is  to  be  stored  in  Mass  Memory  press  "1  Execute." 
The  raw  data  is  stored.   Below  the  tables  will  be 
printed  the  words  "THE  RAW  DATA  IS  STORED  IN  RAW  DAT 
FILE  #?  Along  with  the  file  number  that  was  input  in 
step  12.   If  data  is  not  to  be  stored  press  "0  EXECUTE" 
and  the  words  "THIS  DATA  WAS  NOT  STORED"  will  be  printed. 

20.  The  calculator  will  automatically  CHAIN  with  "TTR  2." 

21.  Temperatures  of  the  nozzle,  labyrinth  and  T   will  be 
printed  in  degrees  Rankine.   The  labyrinth  seal  flow 
rate,  its  Reynolds  number  and  kinetic  energy  factor 
will  be  printed.   Also,  the  nozzle  discharge  coefficient, 
nozzle  flow  rate,  its  Reynolds  number,  and  flow  rate 
will  be  printed. 
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22.  The  display  will  read  "CALC  NEW  W-DOT?  YES=1,  NO=0." 
If  you  elect  to  calculate  a  new  theoretical  flow  rate, 
based  on  a  stator  blockage  factor  of  0.965,  and  employ 
this  value  in^  the  remainder  of  the  program,  press  "1 
EXECUTE."   However,  if  you  use  the  original  flow  rate, 
press  "0  EXECUTE."   Go  to  step  24. 

23.  The  election  for  a  new  calculation  of  flow  rate  was 
chosen;  a  new  flow  rate,  the  stator  flow  function  and 
blockage  factor  will  be  printed  out.   Note  that  the 
flow  function  and  blockage  factor  were  employed  to 
calculate  the  new  flow  rate.   Go  to  step  25. 

24.  The  original  flow  rate  is  used.   The  stator  flow  function 
and  blockage  factor  based  on  the  flow  rate  and  measure- 
ment are  printed. 

25.  The  display  will  read  "INPUT  ALPHA  1  (DEG) ,  OR  0." 
Inputing  any  number  other  than  0  will  fix  the  stator 
exit  flow  angle  in  degrees.  This  method  assumes  VU, 
is  calculated  from  the  stator  torque  and  adjusts  the 
value  of  VA,  to  agree  with  the  elected  stator  exit  angle. 
By  pressing  "0  EXECUTE"  the  stator  exit  angle  is  not 
chosen.   Go  to  step  27. 

26.  By  choosing  to  fix  the  stator  exit  flow  angle,  the 
value  you  selected  will  be  printed.   The  value  of  K, 
which  is  the  ratio  of  VA,  determined  by  the  momentum 
equation  to  VA,  determined  by  continuity,  will  be 
printed.   VAlues  from  1.0  to  1.2  are  within  acceptable 
limits.   Go  to  step  29. 
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27.  The  value  of  the  stator  exit  as  calculated  from 
measurements  will  be  used.   The  display  will  read  "ENTER 
K."   Values  between  1.0  and  1.2  are  within  acceptable 
limits.   Press  " (K)  EXECUTE."   This  value  will  be 
printed  and  the  program  will  CHAIN  with  "TTR  3." 

28.  TTR  3  will  evaluate  Control  Volume  B  and  will  CHAIN 
with  "TTR  4." 

29.  The  display  will  read  "CALC  THEOR.  LOSSES  ?  YES  =  1, 
NO  =  0 . "   This  program  option  permits  you  to  calculate 
both  the  theoretical  rotor  and  stator  loss  coefficients. 
Copies  of  REF.  16  and  18  will  be  needed.   By  pressing 

'1  EXECUTE,"  you  have  elected  to  calculate  the  losses. 
Press  "0  EXECUTE"  if  you  do  not  desire  to  calculate 
the  losses.   Go  to  step  37. 

30.  The  display  will  read  "DEL  ALPHA  =  (SOME  NUMBER) ." 
Note  the  number  and  press  "CONTINUE  EXECUTE." 

31.  The  display  will  read  "BETA  1  =  (SOME  NUMBER)."   Note 
the  number  and  press  "CONTINUE  EXECUTE." 

32.  The  display  will  read  "INPUT  ZETA  PO,  FIG.  15,  1174 
VA-j?"   Turn  to  Fig.  15  in  Ref.  16.   Select  the  value  of 
ZETA  PO  which  corresponds  to  the  values  of  DEL  ALPHA, 
BETA  1;  and  ALPHA  EXIT.   press  ("ZETA  PO)  EXECUTE." 

33.  The  display  will  read  "RE  ROTOR  =  (SOME  NUMBER)." 
Note  this  number  and  press  "CONTINUE  EXECUTE." 

34.  The  display  will  read  "INPUT  K-RE,  FIG.  23,  GA1074VA2." 

Turn  to  Fig.  23  in  Ref.  18  and  select  the  value  of  K_,„ 

Kb 

that  corresponds  to  the  rotor  blade  Reynolds  number 
given  in  step  33.   Press  "(K-RE)  EXECUTE." 
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35.  The  display  will  read  "M2  (IS)  =  (SOME  NUMBER)." 
Note  this  number  and  press  "CONTINUE  EXECUTE." 

36.  The  display  will  read  "INPUT  KM,  FIG.  18,  1174VA1." 
Turn  to  Fig.  18  in  Ref.  16,  and  select  the  value  of  KM 
that  corresponds  to  the  blade  exit  Mach  Number  as 
calculated  from  isentropic  conditions.   Press  "(KM) 
EXECUTE." 

38.  The  program  will  CHAIN  with  "TTR  6."   The  display  will 
read  "ENTER  NEXT  RED  REC.  #."   Input  the  same  record 
number  from  step  12  and  press  "EXECUTE." 

39.  The  display  will  read  "STORE  RED.  DATA  ?  YES  =  1, 

NO  =  0."   Press  "1  EXECUTE"  if  you  desire  reduced  data 
to  be  stored  and  used  for  plotting  and  tabulation. 
Press  "0  EXECUTE"  if  you  do  not  desire  to  store  reduced 
data.   Go  to  step  41. 

40.  You  have  elected  to  store  the  reduced  data  in  the  Mass 
Memory,  and  the  words  "REDUCED  DATA  IS  STORED  IN  FILE/ 
REDDAT/RECORD  (NUMBER  FROM  STEP  12)."  will  be  printed. 

41.  Electing  not  to  store  the  reduced  data  the  words  "THIS 
DATA  WAS  NOT  STORED"  will  be  printed. 

42.  The  display  will  read  " (CONT) (EXEC)  MORE  DATA  -  OR 
(END) ." 

43.  If  more  data  points  are  to  be  reduced,  press  "CONTINUE 
EXECUTE."   The  display  will  read  "TAB  RAW  DATA  ?  YES  =  1 
NO  =  0."   If  you  desire  not  to  tabulate  raw  data,  press 
"0  EXECUTE."   Then  the  display  will  read  "TAB  REDUCED 
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DATA  ?  YES  =  1,  NO  =  0."   If  you  desire  not  to  tabulate 
reduced  data,  press  "0  EXECUTE."   By  inputting  zeros 
for  both  of  these  tabulation  questions  the  program  will 
display  "REDUCE  NEXT  POINT,"  and  will  CHAIN  with  "TTR 
1."   Go  to  step  12.   Note:   Tabulation  is  usually 
elected  after  all  points  are  reduced  for  a  run. 

44.  To  tabulate  the  raw  data  press  "1  EXECUTE."   The  program 
will  GET  "TTR  7."   The  display  will  read  "EXECUTE  RECORD 
#'S:   LOWEST,  HIGHEST."   Input  your  first  and  last 
record  numbers  and  press  "EXECUTE."   The  display  will 
read  "OMIT  4  RECORDS.   ENTER  #'S  OR  0fS."   This  option 
permits  you  to  omit  any  points.   Input  the  record 
numbers  and  press  "EXECUTE." 

45.  The  display  will  read,   read  "REQUEST  DATA  CH  # :   0,1 

OR  10."   Press  "1  EXECUTE"  will  cause  the  4  8  measurements 
to  be  tabulated  under  the  heading:   TTR  INPUT  DATA. 
Afterwards,  the  calculator  will  request  the  next  channel 
number.   Press  "10  EXECUTE"  will  cause  all  other  channels 
to  be  tabulated  and  printed  with  the  same  heading.   Next, 
press  "0  EXECUTE."   The  display  will  read  "TAB  REDUCED 
DATA  ?  YES  =  1,  NO  =  0."   Press  "1  EXECUTE"  and  the 
program  will  GET  "TTR  8."   Go  to  step  46.   Press  "0 
EXECUTE"  and  the  program  will  GET  "TTR  1."   Go  to  step 
12. 

46.  By  electing  to  tabulate  the  reduced  data,  the  program 
will  GET  "TTR  8."   The  display  will  read  "ENTER  LOWEST, 
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HIGHEST  REC.  # . "   Input  the  first  and  last  record 
numbers,  and  press  "EXECUTE." 

47.  The  display  will  read  "ENTER  1ST  AND  LAST  PT.  #." 
Input  the  two  points,  and  press  "EXECUTE.' 

48.  The  display  will  read  "ENTER  RUN  #."   Input  the  run 
number,  and  press  "EXECUTE." 

49.  The  velocities,  Mach  Number,  and  losses  will  be  tabu- 
lated and  printed.  The  tabulation  will  stop  to  allow 
the  page  to  be  removed  from  the  printer.   Then  press 

" CONTINUE  EXECUTE . " 

50.  Tabulation  and  printing  will  resume  until  "DATA  TABULATION 
COMPLETE"  is  printed.   The  display  will  read  "CHAIN 

WITH  TAPE  READER."   There  is  a  3  second  pause ,   and  the 
program  will  GET  "TTR  1."   Go  to  step  12. 
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TABLE  B-2 
DEFINITION  OF  VARIABLES  AND  EQUATIONS 

A  Array 

AO  Area  of  Control  Volume  (Fig.  B-2,  B-3) 

Al  Area  of  Control  Volume 

A3  Alpha  1 

A4  Alpha  2 

A5  Decision  Variables 

A8  Rotor  Shroud  Pressure  Integrated  Over  Area 

A9  Rotor  Shroud  Pressures  Integrated  Over  Area 

BO  EQ  A  (8)  Ref.  5 

Bl  B2  Eq.  A  (16)  Ref.  5 

B2  Beta  1 

B3  Beta  2 

B4  Units  (2*CP*g) 

B5  Rotor  Shroud  Pressures  Integrated  Over  Area 

B6  Rotor  Shroud  Pressures  Integrated  Over  Area 

B7  Rotor  Shroud  Pressures  Integrated  Over  Area 

B8  Del  Alpha 

B9  Mixing  Loss  Coefficient,  Referred  to  Average  Kinetic 
Energy,  Eq.  71,  Ref.  16 

C  Array 

CO  Gas  constant   (ft-lbf/lbm-R) 

CI  Specific  Heat  (BTU/lbm-R) 

C2  Nozzle  Diameter  (in.) 
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2 

C3  Gravitational  Constant   (32.174  lbm-f t/lbf-sec  ) 

C4  Pipe  Diameter   (upstream) 

C5  Ratio  (Dnoz/Dp) 

*   4.      0.16384 

C6  Conversion  factor:   — ■ 

/2.036 

C7  Discharge  Coefficient 

C8  Thermal  Expansion  Coefficient 

C9  Compressibility  Coefficient 

D  Array 

DO  See  Table  B-l 

Dl  See  Table  B-l 

D2  See  Table  B-l 

D3  See  Table  B-l 

D4  See  Table  B-l 

D5  See  Table  B-l 

D6  See  Table  B-l 

D7  See  Table  B-l 

D8  See  Table  B-l 

D9  See  Table  B-l 

E0  Effective  Total-to-Static  Efficiency  Eq.  A  (18)  Ref .  5 

El  Next  record  Number 

E2  Input  Variables 

E3  See  TRble  B-l 

E4  Profile  Loss  Coefficient  of  Turbine  At  Design  Point 

E5  Profile  Loss  Coefficient  (C  )  Eq.  46,  Ref.  16 

E6  Secondary  Flow  Loss  (s  )  Eq.  63,  Ref.  16 

S3 

E7  Mixing  Loss  Coefficient  (<;  )  Eq.  74,  Ref.  16 
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E8       Stator  Loss  Theoretical  (Total  Blading  Losses) 
Eq.  77,  Ref.  16 

E9       Rotor  Loss  Theoretical  (Total  Blading  Losses) 
Eq.  11,    Ref.  16 

FO       Overall  Axial  Force 

Fl       Partial  Sum  of  Overall  Axial  Force 

F2       P  nA  , 
cl  cl 

a  tot 
F4       Stator  Axial  Force 

F5       Closure  Plate  Force 

F6       Ratio:  (Thickness  to  Spacing)  Eq.  73,  Ref.  16 

F7       Mixing  Loss  Coefficient,  Referred  to  Average  Kinetic 
Energy,  Eq.  71,  Ref.  16 

F9       Correction  Factor  for  Blade  Thickness,  Eq.  47,  Ref.  16 

G        Array 

GO        1^-1 

Y 

Gl       Input  Variable 

G2        — ^r 
Y  -  1 

G6  Decision  Variable 

G7  Decision  Variable 

G8  Secondary  Flow  Loss,  Eq.  63,  Ref.  16 

G9  Thickness  to  Spacing,  Eq.  73,  Ref.  16 

H  Pressure  at  Hub 

HO  See  Table  B-l 

HI  Horsepower 

H2  Referred  Horsepower 

H5       Correction  Factor  For  Effects  of  Reynolds  Number 
(KRe> 
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H6       Correction  Factor  For  Blade  Thickness  (K. )  Eq.  47, 
Ref.  16 

H8       Density  Appendix  B,  Ref.  5 

H9       Reynolds  Number  of  Rotor,  Table  B-2,  Ref.  5 

I  First  Point  (Used  in  Tabulation  Program) 

10       Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  DO  in  Table  B-l 

II  Ratio:   Non-Dimensional  Pressures  at  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  Dl  in  Table  B-l 

12  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D2  in  Table  B-l 

13  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D3  in  Table  B-l 

14  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D4  in  Table  B-l 

15  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D5  in  Table  B-l 

16  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D6  in  Table  B-l 

17  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D7  in  Table  B-l 

18  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D8  in  Table  B-l 

19  Ratio:   Non-Dimensional  Pressures  At  Various  Tap 
Locations  In  Stator  (Fig.  5)  -  D9  in  Table  B-l 

J  Mach  #  At  Station  9  In  Stator 

JO  Joules  Constant  (778.16  Ft-lbf/BTU) 

J2  See  Table  B-l 

J4  Input  Variable 

J5  Profile  Loss  Coefficient  Eq.  46,  Ref.  16 

J6  Velocity  at  Station  9  in  Stator 

J7  Non-Dimensional  Stator  Pressure 
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J8  Velocity  at  EX-TIP 

J9  Velocity  at  Station  5  in  Stator 

K  Fourth  Degree  Polynomial  For  K.E  Factor   (Fig.  4) 

KO  Viscosity  (Sutherland  Formula) 

Kl  Kinetic  Energy  Factor   (G.  E.  Report) 

K2  Beta   (G.  E.  Report) 

K3  K  x  B  x  h^   x  PL  (KxK2xk6xQ0) 

K4  Isentropic  Head  Coefficient  Eq.  a  (22) ,  Ref .  5 

K5  Slope  of  K.  E.  Curve  (Newtonian  Method) 

K6  Area  of  Labyrinth 

K7  K.  E.  Factor  (Convergence) 

K8  Second  Degree  Polynomial  K  (R  )  For  Nozzle  (Fig.  5) 

K9  See  Table  B-l 

L  Counter 

LO  X2  a2  +  (Xu2  -  U2)2  Eq.  A  (21),  Ref.  5 

LI  Rotor  Loss  Coefficient  Eq.  A  (21) ,  Ref.  5 

L2       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L3       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L4       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L5       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L6       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L7       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

L8       Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 
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L9  Mach  #  At  Various  Pressure  Tap  Locations  in  Stator 
(Fig.  5) 

MO  Referred  Moment 

Ml  Mach  #  At  Station  1 

M2  Axial  Mach  #  At  Station  1 

M3  Mach  #  At  Station  2 

M4  Axial  Mach  #  At  Station  2 

M5  See  Table  B-l 

M6  See  Table  B-l 

M7  Referred  Stator  Moment 

M8  Mach  #  At  Exit  -  Tip 

M9  Mach  #  At  SS  -  Aft 

NO  See  Table  B-l 

Nl  See  Table  B-l 

N2  See  Table  B-l 

N3  See  Table  B-l 

N4  See  Table  B-l 

N5  See  Table  B-l 

N6  See  Table  B-l 

N7  See  Table  B-l 

N8  See  Table  B-l 

N9  See  Table  B-l 

00  Incidence  Loss  Coefficient 

01  See  Table  B-l 

02  Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

03  Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 
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04 

05 

06 

07 

08 

09 

PO 
PI 

P2 

P3 

P4 

P5 

P6 

P7 

P9 

QO 

Ql 

Q2 

Q3 

Q4 

Q5 

R 

Re 


Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Velocities  at  Various  Pressure  Tap  Locations  In 
Stator  (Fig.  5) 

Reference  Pressure   (29.92  in.  Hg) 

P-HOOD/P-LAB (RATIO) 


o 


X 


ul 


(sin  a,) 


■]/x 


al 


P,/P   Eq.  A  (10)  Ref.  5 
1     ro 

P2 

P1/Pt1 

Pressure  Ratio  (P  /V2) 

See  Table  B-l 

Average  Pressure  At  Rotor  Exit 

See  Table  B-l 

See  Table  B-l 

P  /P 
to'  atm 

See  Table  B-l 

See  Table  B-l 

Decision  Variable 

Array 

Reynolds  Number  (Nozzle) 


(Alpha  Input) 
(K  Input) 
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RO  Referred  RPM  (N//6) 

Rl  Reynolds  Number  (Labyrinth) 

R2  REq.  A  (4),  Ref .  5 

R3  Theoretical  Degree  of  Reaction  Eq.  A  (23) ,  Ref.  5 

R4  Actual  Degree  of  Reaction  Eq.  A  (24)  ,  Ref.  5 

R8  See  Table  B-l 

R9  Ratio  (14 . 696/Patm) 

S  Array 

S  Temperature  Subroutine  Variable 

50  p1 

51  Temperature  Value 

52  Temperature  Value 

53  See  Table  B-l 

54  See  Table  B-l 

55  See  Table  B-l 

56  See  Table  B-l 

57  See  Table  B-l 

58  See  Table  B-l 

59  See  Table  B-l 

T  Pressure  at  Tip 

TO  Reference  Temperature   (518.7  °R) 

T2  T     (Upstream  Temperature  °R) 

T3  /e"  Temperature  Ratio 

T5  T    (Nozzle  Temperature  °R) 

T6  T    (Plenum  Temperature  °R) 

T7  Te  (Fig.  B-4) 


T8 


Tx  (Fig.  B-4) 
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T9  T2is  (Fig.  B-4) 

U  Array 

UO  Mixing  Loss  Coefficient,  Eq.  74,  Ref.  16 

Ul  (DQ  +  Di)/48 

U2  (DQ  +  Di)/48 

U3  Dimension less  Velocity 

U4  Dimensionless  Velocity 

U5  Factor  %,    P.  37,  Ref.  16 

U6  Loss  Due  to  Tip  Clearance  Flow,  Eq.  75,  Ref.  16 

U7  Tip  Clearance  Loss  Coefficient,  Referred  To  Isen- 
tropic  Conditions  Eq.  76,  Ref.  16 

08  Pto/pi 

U9  Stator  Flow  Function 

V  Array 

Tangential  Rotor  Velocity  (V  ■,) 


VO 


VI  Limiting  Velocity  (V  ) 

V2  V  2  Average  Tangential  Velocity  Eq.  A  (11) ,  Ref.  5 

V3  Val 

V4  Va2 

V5  V2 

V6  V± 

V7  Eq.  A  (10) ,  Ref.  5 

V8  Stator  Blockage  Factor 

W  Array 

WO  Total  Flow  Rate  (W  =  W    -  W,  .  ) 

noz    lab 

Wl  Labyrinth  Flow  Rate  Based  on  K.  E.  Factor  (G.  E. 
Report) 
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W2       Labyrinth  Convergence  Flow  Rate 
W3       Flow  Rate  (ASME  Equation) 
W4       Ratio  (WO/g) 
W5       Referred  Flow  Rate 
W6       W^^  (See  Fig.  B-5) 
W2  (See  Fig.  B-5) 


W7 


W8       W*   =  W,  Cos  {B1   -  62) 
W9       W  2  (See  Fig.  B-5) 
X        Array 


XO 


Dimensionless  Velocity  X  ■, 

1  -  Eq.  A  (7) ,  Ref .  5  (K  Input) 


XI  X 

XI  X  ,  (Alpha  Input) 

X2  X  =  X  0  +  X  2 
u2     a 

X3  1/A  PdA 

X4  X  2  Eq.  A  (7) ,  Ref.  5 

X5  X^  Eq.  A  (15)  ,  Ref.  5 

X6  X2 

X7  R/PtQA   (Alpha  Input) 

X8  Y(2K-1)  +  1 

X9  Denominator  of  Eq.  63,  Ref.  16 

Y  Array 


Hi 

YO        /  1  "Cp^-)  y  For  10 

to 


Yl        /  1  -(^— )  y  For  II 

to 


1EP 

Y2  /  1  "Cn1-)  Y        For  I2 


V      to 
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F- 


/  (f-)       Y  For  13 

to 


/ 


Y4  /      1   -    (-£-)       Y  For   14 

Fto 


A7 


P       Hi 
Y5  ^    /    1        CpM       Y  For   15 

/  V  *to 


F 


P       S 
Y6  /    1   -    (=£-)       Y  For   16 

to 


/ 


— p — ET 

I  1   -    C~-)       Y  For  17 

to 


/ 


P        1^ 
Y8  /l  -   (^— )       Y  For  18 

to 


P      Izi1 

Y9  /l   -    (——)        Y  For   19 

to 


/ 


Z  Array 

ZO  CI  Eq.    A    (8a) ,    Ref.    5 

Zl  Stator  Loss  Coefficient 

Z2  C2  Eq.  A  (16) ,  Ref.  5 

Z3  Rotor  Loss  Coefficient  Eq.  A  (21) ,  Ref.  5 

Z6  Non-Dimensional  Pressure  At  Exit  -  Hub 

Z7  Loss  Due  Tip  Clearance,  Eq.  75,  Ref.  16 

Z8  Rotor  Loss  Coefficient   (Other  Factors) 
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FIG.  B-l  DATA  REDUCTION  SCHEMATIC 


TTR  1 

1.) 

RAW  DATA  STORAGE 

2.) 

RAW  DATA  PRINTED 

3.) 

CHANNEL  ASSIGNMENT 

4.) 

TEMPERATURE  SUBROUTINE 

\t 

TTR  2 

1.) 

CALCULATE  MASS  FLOW  RATE 

2.) 

EVALUATION  OF  CONTROL 
VOLUME  A  (  FIG.  B-2) 

* 

r 

TTR  3 

1.) 

EVALUATION  OF  CONTROL 
VOLUME  B  AND  TEMPERATURES 
(  FIGS.  B-3,  B-4) 

* 

' 

TTR  4 

U) 

CALCULATE  STATOR  EXIT 
VELOCITIES  (FIG.  B-5) 

2.) 

CALCULATE  THEORETICAL 
LOSSES 

i 
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FIG.    B-l  DATA  REDUCTION  SCHEMATIC 


J, 

TTR  5 

1.) 

PRINT  REDUCED  DATA 

\f 

TTR  6 

1.) 

STORE  REDUCED  DATA 

> 

f                                     - 

1 

TTR  7 

1.) 

TABULATION  OF  RAW  DATA 

v 

/ 

TTR  8 

I.) 

TABULATION  OF  REDUCED  DATi 
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FIGURE  B-2       CONTROL  VOLUME   a 


STATOR 


FIGURE  B-3   CONTROL  VOLUME  b 
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FIGURE  B-4  THERMODYNAMIC  PROCESS  OF  FLUID 
IN  AN  AXIAL  TURBINE  STAGE 
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